ONDURWNPE 


. LOi noi dau - Gido trinh thuc hanh Cad 
. Thao tac trong cUta s6 1énh cUla matlab 
. Ham va Script files 

. symbolic va Simulink 

. M6 hinh hé th6ng vién thong 


. Tao giao dién trong Mathlab 
. Thiét ké - m6 phong mach va vé mach in 


LOi ni dau - Gido trinh thuc hanh Cad 
LOi ndéi dau 


Cing v6i sy phat trién nhanh chong cUa may tinh, CAD (Computer-Aided 
Desgin) duc xay dung ngay cang hoan thién va Ung dung trong hau hét 
cac linh vc khoa hoc ky thuat. D6i vGi chuyén nganh Dién ttt, nhiéu 
phan mém CAD cho phép thiét ké mach, m6 phong va vé mach in mot 
cach nhanh chong va hiéu qua nhu OrCAD/Pspice, Multisim (Electronics 
WorkBench), MicroSim, ExpeditionPCB, ... 


Tuy nhién, day 1a cc phan mém dong goi chi dudc Ung dung trong 
chuyén m6n hep 1a Dién tu, n6é khéng cho phép lap trinh m6 phOng cac hé 
thong déng (Dynamic systems) bat ky khac. Vi vay, chUOng trinh dao tao 
mon CAD cho sinh vién Dién ttf chuyén nganh Vién thong va Tu dOng 
héa da huGng tdi phan mém Matlab. Day 1a mOt ng6n ngf lap trinh cap 
cao dang ngu6n m@, no hé tro rat nhiéu thu vién chUc nang chuyén biét 
tU’ Todn hoc, Kinh té, Logic m0, Truyén théng, Diéu khién tu dOng, ... 
dén diéu khién phan cUing cho cdc thiét bi. Dong thdi, né cho phép ngudi 
sur dung b6 sung cac cong cU tU tao lam phong phi thém kha nang phan 
tich, thiét ké va m6 phOng cdc hé thong d6ng lién tuc va rdi rac, tuyén 
tinh va phi tuyén bat ky. V6i nhting uu diém ndi bac cUa minh, Matlab 
da duc nhiéu trudng Dai hoc hang dau trén thé gidi 4p dung va giang 
day. 


Tham voOng thi nhiéu nhung trong pham vi 30 tiét thUc hanh, ching ta chi 
c6 thé dé cap dén nhttng van dé co ban nhat. Hy vong tU do sinh vién 
tu nghién ctu, hoc hoi dé cé thé thiét ké va m6 phOng cdc m6 hinh Hé 
thong Vién thong hay cdc Hé diéu khién ty d6ng. D6ng thdi cing nam 
sO ludc vé cdc phan mém m6 phong mach va vé mach in nhu Multisim, 
OrGad... 


Gido trinh nay g6m 7 bai, m6i bai 5 tiét. Sinh vién chon 6 bai dé thuc tap: 
1. Thao tac trong cUfa sO 1énh cla Matlab. 


2. Ham va Script files. 


3. Symbolic va Simulink. 

4. M6 hinh hé théng Vién théng.(sinh vién chuyén nganh VT) 

5. M6 hinh hé thong Diéu khién tu dOng.(sinh vién chuyén nganh DKTD) 
6. Tao giao dién trong Matlab. 


7. Thiét ké - M6 phOng va vé mach in(cac sinh vién khong chon bai 4 
hoac 5) 


Mac di da hét sttc co gang, song do trinh dO han ché cua nguoi viét ma 
nhiéu van dé chac chua du@c trinh bay tot cting nhu chua bo cUc hgp ly. 
Xin chan thanh cam On mdi y kién déng gop cua sinh vién va cdc ban 
d6ng nghiép. 

TcAD, thang 11 nam 2003 

Nguyén Chi Ng6on 

Dia chi lién hé: 

BO mon Vién thong va Tu d6ng hda 

Khoa Céng nghé Théng tin, Dai hoc Can tho 

01 Ly Tu Trong, Tp. Can tho, tinh Can tho 

Tel: (71) 831301Fax: (71) 830841 

Email: ncngon@ctu.edu.vn 

URL: http://www.cit.ctu.edu.vn/department/ac/ncngon.html 

Muc luc 

LOi ndi daul 


Muc luc3 


Thao tac trong cUta sO 1€nh cUa matlab 
Day 1a gido trinh vé thao tac trong cUfa sO 1@nh cua Matlab 


MUc tiéu 
Bai thi nghiém nay gitp sinh vién lam quen nhanh vGi Matlab 6.5 bang cac 


thao tac don gian trén ma tran, vecto, biéu thtc todn hac, cac 1énh d6 
hoa, ..., thuc hi€n ngay trén cUfa s6 1énh (command window) cUa Matlab. 


Tham khAo 


[1].Nguyén Htta Duy Khang, Bai gidng mén CAD, BO mon Vién thong & 
Tu dOng hda, khoa Cong nghé thong tin, Dai hoc Can tho, 2001. 


[2].The Mathworks Inc., Matlab Notebook User’s Guide, 2003. 


[3].Nguyén Hoai Son - D6 Thanh Viét - Bui Xuan Lam, U’ng dung 
MATLAB trong tinh todn ky thuat, Tap 1, NXB DHQG Tp. HCM, 2000 


[4]. Nguyén Httu Tinh - Lé Tan Hing - Pham Thi Ngoc Yén - Nguyén 
Thi Lan Huong, Co sO Matlab & Ung dung, NXB KH va ky thuat, 1999. 


[5].http://www.facstaff.bucknell. edu/maneval/help211/exercises.html 


[6]. http://www. glue.umd.edu/~nsw/ench250/matlab.htm 


ThUc hanh 


TU cUfa sO lénh cla Matlab, sinh vién lan luot thuc hién cac thao tac sau: 


Ma tran 


Dé tao ma tran trong Matlab ta chi can liét céc phan tU’ cUa ma tran trong 
cap dau ngoac vuong ([...]). Cac phan tUf trén cting hang duc phan biét 
bdi dau phay (,) hodc khoang trang (space). Cac hang cUa ma tran, phan 


cach nhau boi dau cham phay (;). Vi du, nhdap ma tran A co 4 hang, 4 ct 
nhu sau: 


>> A=[163213;510118;96712;415141] 
>> size(A) 


Dé truy xudt dén ting phan tl cUa ma tran ta ding chi s6 phan tr tuong 
Ung. Vi du, phan tt’ hang thf 2, cét tht’ 3 cua A 1a A(2,3). 


>> A(2,3) 
Cho ma tran A=[2 4 1; 6 7 2; 35 9], sinh vién ding cdc lénh can thiét dé: 


e Lay dong dau tién cUa ma tran A. 

¢ Tao ma tran B bang 2 dong cudi cing cla A. 

¢ Tinh tOng cdc phan tUr trén cdc ct cUa A. (gQi ¥: tinh tong cdc phan 
tU trén cOt 1: sum(A(:,1))). 

¢ Tinh tOng cdc phan tll trén cdc dong cUa A. 


Cho ma tran A=[2 79 7;3156; 8125], sinh vién giai thich két qua 
cUa cac 1énh sau: 


A 
. AC;[1 4]) 

- A([2 3],[3 1]) 

. reshape(A,2,6) 

A(:) 

. LA A(end,:)] 

. A(1:3,:) 

.[A; A(1:2,:)] 

sum(A) 

10. sum(A') 

11. [| A; sum(A) | [ sum(A,2) ; sum(A(:)) J J 


WOANAUARAWNHE 


Giai hé phuong Ax=b, v6i: A= 2 5 3 vwab= 1 .GOiY¥: 


x=A\b. 


VectO 


Vecto thUc chat cting 1a ma tran c6 kich thuGc (n x 1) hay (1 x n), nén ta 
c6 thé tao ra vectO nhu cach tao ra ma tran. Ngoai ra, cd thé ding m6t s6 
cach sau: 


>>x=0:0.1:1 
>>y=linspace(1, 10, 20) % vecto 20 phan tu cach deu nhau tu 1 den 10 
>>z=rand(10,1) 

Cho vectO x = [315 79 2 6], gidi thich két qua cUa cac lénh sau: 


e x(3) 

e x(1:7) 

e x(1:end) 

e x(1:end-1) 

e x(6:-2:1) 

e x({(1 621 1)) 
° sum(x) 


Sinh vién hay tao m6t vecto x cd 100 phan ttt, sao cho: 


_ (—1)"1 
in = “n-I) * 
Gdi y: Tao vecto n cé 100 phan tu tl 1 dén 100, ding toén tU dau cham 
(.) dé xac dinh x. 


Cac da thUtc 

Cac da thttc trong Matlab dugc m6 ta bang cac vecto hang v6i cdc phan 
tU cUa vecto chinh 1a cdc hé s6 cUa da thc, xép theo tht tu's6 mii gidm 
dan. Vi du, da thc m = s4-s3+4s2-5s-1 duoc biéu dién 1a: 

>>m=[1 -145 -1] 


Dé xac dinh gia tri cUla da thUtc, ta dting 1énh polyval. Vi du, xac dinh gid 
tri cUa da thtrc tai diém s=2: 


>>polyval(m,2) 
Dé xac dinh nghiém cua da thc, ta ding 1énh roots. Vi du: 
>>roots(m) 

Cho phuong trinh x2-4x+5=0, giai phuoOng trinh theo 2 cach, cach 1 — 
tinh delta theo phu‘ong phap cO dién, cach 2 — dung ham roots, hay so sanh 
két qua. 

¢ Cachl: 
>>a=1; 
>>b=-4; 
PaeH=o 
>>x1=(-b+sqrt(bA2-4*a*c))/(2*a) 
>>x2=(-b-sqrt(bA2-4*a*c))/(2*a) 
e Cach 2: 
>>m=[a b c]; 


>>x=roots(m) 


Hay thay d6i cdc gid tri khac nhau cUa a, b va c tuOng Ung trong 2 cach 
giai trén. So sanh két qua va nhan xét. 


Giai phuOng trinh x3- 2x2+4x+5=0. Kiém ching két qua thu dugc bang 
ham polyval. Sinh vién cé nhan xét gi vé két qua kiém chUng. 


Lap lai cau cho phuoOng trinh x7-2=0. 
Sinh vién tht ding ham poly dé tao da thtfc tl’ céc nghiém cho tru6c. 


>>help poly 


D6 hoa 


Matlab h6 tro ché dO dé hoa rat manh, bao g6m d6 hoa 2D va 3D, voi 
cdc truc toa d6 tuyén tinh va phi tuyén bat ky. 


D6 hQa 2 D 
DO hoa 2D chu yéu dya trén 1énh plot. Dé dugc gitip dG, ta go: 
>>help plot 


Vé dé thi ham s6 y1=sinx.cos2x va ham s6 y2=sinx2 trong [0-2 ], trén 
cung hé truc toa dd, ta lan uot thUc hién nhu sau: 


>>x=0:0.01:2*pi; 

>>y1=sin(x).*cos(2*x);%enhan tuong tung tung phan tu 
>>plot(x,y1) 

>>grid on%hien thi luoi 


Sau khi thu duc d6 thi ham y1, dé vé y2 trén cing d6 thi, ta thc hién: 


>>hold on%giu hinh, mac nhien la hold off 
>>y2=sin(x.42);%luy thua tung phan tu 
>>plot(x,y2,’k’)%duong ve co mau den 

>>axis([O 4*pi —1.25 1.25])%dinh lai toa do hien thi 

Ta co thé dat nhan cho cac truc cling nhu tiéu dé cho dé thi: 
>>xlabel(‘Time’) 

>>ylabel(‘Amplitude’) 

>>title(‘y1=sinx.cos2x and y2=sin(x/2)’) 


>>legend(‘sinx.cos2x’,’sinx/2’) 


[missing_resource: graphics1.wmf] 
Hinh 1.1 — Biéu dién d6 thi cdc ham s6 trén cling hé truc toa dé 


Matlab h6 tro rat nhiéu thudc tinh d6 hoa, dé cé thé kiém soat cac 
thu6c tinh nay ta can ding dén thé d6 hoa. Vi du: 


>>close all 

>>x=[0 1 2 3]; 
>>y=[0 4 1 5]; 
>>h=plot(x,y) 


h chinh 1a thé d6 hoa cUa ham plot, dé thay cdc thu6c tinh d6 hoa, ta 
ding 1@nh: 


>>set(h) 


Bay gid ta thr dat m6t sO thudc tinh dé hoa cho h. 
>>set(h,’Color’, ‘r’)%dat lai mau do 
>>set(h,'LineWidth',6)%dat do rong duong 


>>set(h,'"Marker','v',,MarkerSize',6) 


[missing_resource: graphics2.wmf] 

Hinh 1.2 — Thay d6i thuéc tinh duOng biéu dién 

Dt liéu ding dé vé dude dat trong ‘XData’ va ‘YData’. 

>> set(h,'XData',[0 1 1 3]) 

>>set(h,'Y Data’,[0 3 5 1]) 

TU két qua nay, ta nhan thay néu d(r liéu trong ‘XData’ va ‘YData’ bién 


thién theo thdi gian, thi trén clfa s6 d6 hoa (figure) ta sé thdy hinh anh 
sinh déng kiéu ‘animation’. 


[missing_resource: graphics3.wmf] 
Hinh 1.3 - Thay d6i dt? liéu trong clfa s6 d6 hoa 


TuOng tl nhu ham plot, sinh vién thu ding cdc ham semilogx, semilogy 
va loglog cho truOng hgp trUuc toa dO phi tuyén. 


Ngoai cac 1énh biéu dién dung cong trong toa dO Descartes, Matlab 
cing h6 tro viéc vé d6 thi ham s6 trong hé toa dO cU’c bang ham polar. 


>>theta=0:0.05:2*pi; 


>>r=sin(5*theta); 


>>polar(theta,r) 


[missing_resource: graphics4.wmf] 


Hinh 1.4 — Biéu dién d6 thi ham s6 trong hé toa dO cuc 


D6 hga 3 D 

Matlab cung cap nhiéu ham vé do thi 3D, chang han: plot3 - ding dé vé 

cdc dung trong khong gian 3 chiéu; mesh va surf - ding dé vé vat thé 

3D (g6 help mesh va help surf dé biét thém cac ham 3D cé lién quan). 
Vé d6 thi 3D bang ham plot3: 

>>t=0:pi/50:10* pi; 

>>x=sin(t); 

>>y=cos(t); 

S376 

>>subplot(121), plot3(x,y,z)%ve tren o thu nhat 

>>grid on 

>>subplot(122), plot3(x,y,t.\2)%ve tren o thu hai 


>> grid on 


[missing_resource: graphics5.wmf] 


Hinh 1.5 — Vé d6 thi 3D bang ham plot3 


Vé mat paraboloid z=x2+y2 trong khéng gian 3 chiéu: 
>>close all 
>>t=-5:0.1:5; 
>> [x,y]=meshgrid(t);%dinh luoi ve 
>>Z=xA2t+yA2; 
>> subplot(2,2,1), mesh(z)%ve mat luoi 3D 
>> title(‘mesh(z)') 
>> subplot(2,2,2), meshc(z)%giong mesh nhung co them duong vien 
>> title(‘meshc(z)') 
>> subplot(2,2,3), meshz(z)%co them luoi tren mat x,y 
>> title(‘meshz(z)') 
>> subplot(2,2,4), waterfall(z)%chi ve luoi theo 1 huong 
>> title(‘waterfall(z)') 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


Hinh 1.6 - Vé mat paraboloid 


P, iS j 2 2 = 2 coms 
Vé mat z = ae trong khéng gian 3 chiéu: 
m+y 
>>x=-8:0.5:8; 
ee ies 


>>[x,y]=meshgrid(x,y); 


>>r=sqrt(x.A2+y.42); 


>>z=sin(r)./T; 


>>surf(x,y,Z) 


[missing_resource: graphics6.wmf] 
Hinh 1.7 — MOt biéu dién d6 thi 3D khac 


Sinh vién tht’ vé mat tru z = ./24 + y? bang ham mesh va ham surf. 


TV chOn 

Giai hé phuoOng trinh sau: 
2x1+4x2+6x3 — 2x4=0 
X1+2x2+ x3 + 2x4=1 
2X2+4x3 + 2x4=2 
3x1- x2+ 10x4=10. 


ChUing to rang (A+B)C=AC+BG, vi: 


10 -2 3 41 3 4 
A= 20 4 ,B= -10 2 vaC= | P ‘ 
3. 66 0 5 


Sinh vién thu vé hinh sau (Hinh 1.8): 
**EKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 1.8) 


TU ctra s6 1é€nh cUa Matlab, nhap: demos. Chon MATLAB. Graphics 
3D-plots. Chay chuOng trinh demo nay. 


Ham va Script files 

Day la giao trinh vé ham va script files 

Muc tieu 

Bai thi nghiém nay gitip sinh vién ti€p can vi ky thuat lap trinh trong 
Matlab thong qua viéc su dung cac ham cé san, xay dung cdc ham mdi 
dua trén cdc cau tric if - for - while va tim hiéu cach xay dung script files, 
lam cO s6 cho viéc lap trinh m6 phOng cac hé th6ng Vién thong va TU 
dOng 6 cdc bai sau. 


Tham khao 


[1].Nguyén Htta Duy Khang, Bai gidng mon CAD, BO mon Vién thong & 
Tu dOng hda, khoa Cong nghé thong tin, Dai hoc Can tho, 2001. 


[2].The Mathworks Inc., Matlab Notebook User’s Guide, 2003. 


[3].Nguy€n Cong Dinh, Phan tich va t6ng hop cac hé thOng diéu khién 
bang may tinh, NXB Khoa hoc va Ky thuat, 2002. 


[4]. Nguyén Httu Tinh - Lé Tan Hing - Pham Thi Ngoc Yén - Nguyén 
Thi Lan Huong, Co sO Matlab & Ung dung, NXB KH va ky thuat, 1999. 


[5].http://www.mines.utah.edu/gg_computer_seminar/matlab/ 


[6]. http://www. glue.umd.edu/~nsw/ench250/matlab.htm 


ThUc hanh 


Ham va Script files trong Matlab déu duc quan ly dui dang cac tap tin 
cé phan m6 rOng .m, thuOng duc soan thao bdi Matlab Editor. Khdi 
d6ng Matlab Editor bang m6t trong cdc cach sau: 


e Nhdap chudt vao biéu tugng trén menu bar cUa cUfa s6 1€nh. 
e File New M-file 


e Nhdp vao biéu tuong Matlab Editor trén desktop (néu cd). 


Ching ta lan luot khao sat qui cach xy dung cdc ham va script file. Sinh 
vién doc va thUc hanh lan luot theo cac phan sau. 


Script files 

Tap hgp cac dong 1énh cla Matlab dudc sap xép theo mét cau tric nao 
d6 va lutu thanh file cé phan md rOng *.m dugc goi la script file (file kich 
ban, file chUOng trinh). Ta cé thé chay file nay tu cra s6 1énh gidng hét 
nhu cac lénh cUa Matlab. Cau tric cla m6t script file nhu sau: 

- %Phan viét sau dau ‘%’ 6 day ding cho 1énh help 

%Thong thuOng phan nay m6 ta chUc nang, céch str dung, 

%vi dU minh hoa hay nhttng luu y dac biét ma tac gid mong muOn trod 


% gitip cho nguOi sU dung. 


[global ténbién1, ténbién2,... 1% Khai bdo bién toan cuc % (néu C6) 
<cac cau lénh>% phan trinh bay cau 1énh 


Hay khdi d6ng Matlab Editor va tao m6t script file cé tén bai21.m, v6i 
ndi dung nhu sau: 


% Doan script file nay hien thi loi chao trong 2s. Sau do 
% hien thi logo cua matlab mot cach sinh dong roi thoat 


close all 


figure(‘Color'’,[0 0 O],... 
'Name','Welcome to Matlab Experiments’... 
'NumberTitle’,'off',... 


'MenuBar','none’); 


text(‘String’,'Welcome to MATLAB',... 
‘Color;25 25:25 | ,: 

'Position',[0.01 .501],... 
'Fontsize',32.... 

'FontAngle'’,'italic'); 

text( 'String','Welcome to MATLAB"... 
COlor wes: 

'Position',[O .5],... 

'Fontsize',32.... 

'FontAngle'’,'italic'); 

axis off; 


pause(2); % dung trong 2 giay 


logospin 


% ---- Thoat - xoa cac bien trong workspace va dong cua so lai ----- 


clear 


close 


% ket thuc script file 


Sau khi luu file nay, tt’ cUfa s6 16nh cUa Matlab, sinh vién hay nhap: 


>>help bai21 
Dé thi hanh script file va soan, hay nhap: 


>>bai21 


Soan thao script file c6é tén bai22.m dé lap biéu d6 tOng s6 sinh vién 
Dién TU - Tin hoc t6t nghi€p tai khoa Céng nghé Thong tin, Dai hoc 
C&n tho tt nam 1996 dén 2001 véi dtr liéu nhu sau: 


Nam 


1996 


1997 


1998 


1999 


2000 


2001 


KY su Dién tur 
38 
33 
36 
31 
60 


70 


Ky su’ Tin hoc 
48 

94 

120 

a2 

110 


131 


N6i dung cUa script file nhu sau: 

% 

% Lap bieu do tong so sinh vien Dien tu — Tin hoc tot nghiep tai 

% Khoa Cong nghe Thong tin, Dai hoc Can tho, tu 1996-2001 

% 

svdt=[38 33 36 31 60 70];% sinh vien Dien tu 

svth=[48 54 120 92 110 131];%sinh vien Tin hoc 

nam=1996:2001; 

subplot(211), bar(nam,svdt); 

title(‘Sinh vien Dien tu tot nghiep tu 1996-2001'); 

subplot(212), bar(nam,svth); 

title('Sinh vien Tin hoc tot nghiep tu 1996-2001’); 

Colormap(cool); 

% 

Thi hanh file nay tu cUta sO 1énh. Sinh vién c6é thé ding 1énh help bar dé 
biét cac thong sO cUa ham nay va hiéu chinh (tty thich) m6t sO chUc nang 
hién thi cua ham bar trong file bai22.m. 


e Sinh vién hay tao m6t script file c6 tén bai23.m dé vé cac dé thi trén 
hinh 2.1. 


SU dung cac ham xay dng san 


Matlab h6 tro mOt thu vién ham rat phong phu, xay dUng trén cdc gidi 
thuat nhanh va c6 d6 chinh xdc cao. Ngoai cdc ham co ban cUa Matlab, 
tap hgp cac ham ding dé giai quyét mOt Ung dung chuyén biét nao dé 
goi la Toolbox, vi du: XU ly s6 tin hiéu (Digital Signal Processing), Diéu 
khién tu d6ng (Control), Mang Noron nhan tao (Neural networks), ... 
Sinh vién xem lai gido trinh dé biét thém vé cdc ham. Ngoai ra, c6 thé 
ding 1€nh help dé biét chUc nang cUa toolbox va ham cting nhu cach thtfc 
su dung chting. 

help <ten toolbox>% chuc nang toolbox 

>>help control % liet ke ham cua control toolbox 

help <ten ham>% chuc nang ham 


>>help plot% chuc nang ham plot 


Ta c6 thé tim kiém cdc ham lién quan bang cach cung cap cho ham 
lookfor cua Matlab m6t tU khéa: 


lookfor <tu khoa tim kiem> 

>>lookfor filter % tim cdc ham lién quan dén mach loc 
***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 2.1 — dung cho cau III.1.3) 


Ham [Y I]=max(X) cho biét phan ttr 16n nhat cUa vecto (mang) X vGi 
chi s6 tuOng ttmg I. 


>>help max 
>> x=[159 764]; 


>> [y,i]=max(x) 


2) 

Nghia la phan tU tht? i=3 cUa vecto x cé gié tri 16n nhat, y=9. 
Ham Y=exp(X), tinh eX, két qua tra vé cho Y. 

>>help exp 

>>x=0:0.05:1; 

>>y=exp(x); 

>>stem(x,y)% giong ham plot nhung ve cho tin hieu roi rac 


Gia sU ta muOn diéu ché bién dé sdéng mang dudc truyén hai bang canh 
(Double sideband transmission carrier Amplitude Modulation) v6i: 


¢ Tin hiéu diéu ché X=0.2*sin(2*pi*t) vGi t:0:0.01:10 
e Tans6 séng mang fc=620KHz 


Sinh vién hay tim ham thich hop va vé tin hiéu thu duc sau khi diéu ché, 
co dang nhu hinh 2.2: 


[missing_resource: graphics1.wmf] 


Hinh 2.2 — Tin hiéu thu duoc sau khi diéu ché bién d6 


Xay dung ham 


Viéc xay dung ham ciing duc thuc hién tuong tu nhu script file. Tuy 
nhién, d6i vGi ham ta can quan tam dén cdc tham sO truyén cho ham va 
cdc két qua tra vé sau khi thc hién. C6 3 diém can luu y: 
e Tén ham phai dugc dat tring vGi tén file luu tr. 
e Phai c6 tu khda function 6 dong dau tién. 
¢ Trong mt ham cé thé xay dung nhiéu ham con (diéu nay khéng cé 
trong script file). Két thtic ham con phai c6 tU khéa end (diéu nay 
kh6ng can trong ham ‘cha’). 
Qui cdch xay dung ham duc m6 ta nhu sau: 


function [out1,out2,...]=tenham(in1,in2,...) 


[global <ténbién1, ténbién2, ...>]%khai bao bién toan cUc (néu cé) 
<Céc cau lénh thurc hién ham> 
outl=két qua1%két qua tra vé cua ham 


out2=két qua2 


% Cac ham con (néu cé) 
[ function [subout1,subout2,...]=tenhamcon(subin1,subin2,...) 
<Céc cau lénh cua ham con> 


end | %tu khéa end khong can doi voi Matlab Version 6.x 


Xay dung ham gptb2 dé giai phuong trinh bac hai ax2+bx+c=0. N6i 
dung ham nhu sau: 


function [x1,x2]=gptb2(a,b,c) 

% Giai phuong trinh bac hai ax42+bx+c=0 

% [x1,x2]=gptb2(a,b,c) 

% Trong do: x1,x2 nghiem thuc hoac phuc 

% a,b,c la 3 he so cua phuong trinh 

% 

% Vi du: [x1,x2]=gptb2(1,-3,2) 

% 

% Copyright 2003 Nguyen Chi Ngon TcAD - CIT - Cantho University 
% Email: ncngon@cit.ctu.edu.vn 

if nargin<3 

error('Vui long nhap du 3 he so cua phuong trinh’) 
elseif a==0 

x1=-c/b; 

x2=[]; 

else 

D = bA 2 - 4*a*c; 


x1 = (-b+sqrt(D))/(2*a); 


x2 = (-b-sqrt(D))/(2*a); 
end 
Sau khi luu file nay vGi tén gptb2.m, sinh vién th kiém chttng két qua: 
>>help gptb2 
>>[x1,x2]=gptb2(1,6,-7) 
>>[x1,x2]=gptb2(2,7,14) 
>>[x1,x2]=gptb2(0,4,3) 
>>[x1,x2]=gptb2(1,6) 
Cho biét y nghia cla tlr khda nargin? 
Sinh vién hay viét lai ham nay dé két qua chi tra vé nghiém s6 thuc. 
Xay dung ham vdcongdb(a,m,method) dé vé mét s0 dudng cong dac 
biét trong hé toa dO cc, vGi a 1a ban kinh va m 1a s6 dung cong vé trén 
cling truc toa dO. Trung hgp nay ham khéng tra vé gid tri nén ta khong 
can bién ng6 ra. 
Tuy theo gid tri cua tham s6 ‘method’ ma ta vé d6 thi tuong Ung: 


Néu method = ’Becnulli’: V6 duOng Lemniscat Becnulli: 
r = av/| 2c0s20 | 


in30 
1— SS | 


Néu method = ’Astroit’: V6 duOng Astroit: r = a4/ | 
Néu method = ‘Xoanoc’: Vé duGng xoan 6c: r = cos# + 1 
NOi dung ham nhu sau: 


function vdcongdb(a,m,method) 


% Ve duong cong trong toa do cuc: vdcongdb(a,m,method) 
% method = 'Becnulli' - Ve duong Lemniscat Becnulli: 

% r=a*sqrt(abs(2*cos(2*theta))) 

% 'Astroit' - Ve duong Astroit: 

% r=a*sqrt(abs(1-sin(3*theta)/4)) 

% 'Xoanoc' - Ve duong xoan oc: 

% r=a*cos(theta)+1 

% Voi: a-ban kinh; m-so duong cong ve tren cung he truc 
% Vi du: vdcongdb(0.5, 4, "Becnulli') 

% 

% Copyright 2003 Nguyen Chi Ngon, TcAD - CIT - CTU 
% Email: ncngon@cit.ctu.edu.vn 

if nargin<3 

error('Vui long nhap du 3 thong so cua ham’) 

else 

theta=0:0.01:2*pi; method=upper(method); 

switch method 

case 'BECNULLI' 

r=a*sqrt(abs(2*cos(2*theta))); 


case 'ASTROIT' 


r=a*sqrt(abs(1-sin(3*theta)/4)); 

case 'XOANOC' 

r=a*cos(theta)+1; 

otherwise 

error('Chon: "Becnuli", "Aristod" hoac "Xoanoc") 
end % end of switch 

% ve do thi 

close all; figure(‘Color','w'); 

for k=1:m 

hold on 

rl=r*k; 

mau=[rand(1,1) rand(1,1) rand(1,1)]; 
h=polar(theta,r1); 
set(h,'color',mau,'LineWidth’',2); 

axis equal; 

end % end of for 

hold off; 

axis off 

end % end of if 


Sinh vién hay kiém chting lai hoat dOng cUa ham, vi du: 


>>help vdcongdb 
>>vdcongdb(1,5,’ Becnulli’) 
>>vdcongdb(1,5,’ Astroit’) 
>>vdcongdb(1,5,’ Xoanoc’) 
>> vdcongdb(1,5,’saikieu’) 
>> vdcongdb(5, ’becnulli’) 


PP sas 


Xay dung ham dudoan() dé du dodn két qua sau mOi lan tung m6t xtic 
xac dOng nhat, 6 mat. NOi dung ham nhu sau: 


function dudoan() 

% Du doan ket qua sau moi lan tung ngau nhien mot xuc xac 6 mat 

% Chuong trinh lap lai cho den khi nguoi su dung khong doan tiep 

% 

% Copyright 2003, Nguyen Chi Ngon TcAD - CIT, Cantho University 
% Email: ncngon@ctu.edu.vn 

tiep = 'y'; sai=0; dung=0; 

disp(‘Chao mung ban den voi Casino nay!") 

while(lower(tiep)=='y’) 

doan=input('Moi ban du doan ket qua (1-6):’); 


kqua=tungxx; 


if (doan ~= kqua) 

disp('Xin loi, ban da doan sai!') 

Sai=sait1; 

else 

disp('Xin chuc mung!’) 

dung=dung+1; 

end 

tiep=input('Ban muon choi tiep("y"/"n"):'); 
end 


disp(['Dung ' num2str(dung) ' trong tong so ' num2str(sait+dung) ' lan 
doan']) 


% subfunction -------------- 
function mat = tungxx() 
mat=floor(6*rand(1,1))+1; 
% end 


Sinh vién thUr sr dung ham nay va cho biét su khac nhau giCfa script file 
va ham khong c6 tham s6 vao. 


TV chOn 
Viét chUOng trinh tim nghiém thUc cUa phuOng trinh bac 2. 


Viét chuong trinh in tam gidc Pascal n dong trong man hinh d6 hoa v6in 
dudc nhap tU ban phim. 


symbolic va Simulink 
Day 1a gido trinh vé symbolic va simulink 


MUc tiéu 
Hiéu mt cach co ban vé hai cong cu manh va httu hiéu cla Matlab, dé 
1a Symbolic va Simulink, dé tu d6 sinh vién cé thé tl minh phat huy cdc 


chtfc nang cao hOn cua hai cong cU nay trong tinh todn va m6 phong hé 
thong. 


Tham khAo 
[1].The Mathworks Inc., Matlab Notebook User’s Guide, 2003. 


[2].Nguyén Hoai Son - D6 Thanh Viét - Bui Xuan Lam, U’ng dung 
MATLAB trong tinh todn ky thuat, Tap 1, NXB DHQG Tp. HCM, 2000 


[3]. Nguyén Chi Ng6n, Bai thi nghi€m KY thudt m6 phOng trong Diéu 
khién tu d6ng, BO mén Vién Théng & Tu dOng hoa, khoa Céng nghé 
thong tin, Dai hoc Can tho, 2002. 


[4].Nguyén Cong Dinh, Phan tich va t6ng hop cac hé thOng diéu khién 
bang may tinh, NXB Khoa hoc va Ky thuat, 2002. 


[5].http://www-h.eng.cam.ac.uk/help/tpl/programs/Matlab/S ymbolic.html 


ThUc hanh 
Symbolic va Simulink déu chia thu vién chUc nang rat phong phi, bai thi 


nghiém nay chi c6 thé dé cap dén nhting gi goi 1a co ban nhat. TU do, 
sinh vién c6 thé tu minh nghién cUtu va phat trién tiép. 


Symbolic 


Phién ban m6i nhat cua Symbolic toolbox dudc Mathworks gidi thiéu 
trong Matlab 6.5 vao thang 6-2003. Do 1a m6t thu vién toan hoc kiéu ky 
tu, dudc phat trién tur symbolic Maple cUa truOng Dai hoc Waterloo, 
Canada. Dé cé cai nhin t6ng quat vé cdc chic nang cUla Symbolic, sinh 


vien hay go: 


>>help symbolic 


MOt sO ham théng dung cUa Symbolic: 


Tén 
ham 
diff 
int 
taylor 


det 


numden 


subs 


dsolve 


Chu nang 
Dao ham 
Tich phan 


Khai trién Taylor 


Dinh thutc cua ma 
tran 

TU va mau cua 
phan so 


Thay bién sym 
bang tri sO 


Giai phuong trinh 
vi phan 


Tén 
ham 
fourier 
ifourier 


laplace 


laplace 


ezplot 


ezpolar 


ezmesh 


Chu nang 


Bién d6i Fourier 


Bién d6i Fourier 
nguodc 


Bién d6i Laplace 


Bién d6i Laplace 
nguodc 


Vé ham, plot 


Vé ham, toa d6 
cUc_ polar 

vé mat ludi 
mesh 


solve Giai phuOng trinh ezsurf Vé mat surf 
dai sO 


Dé bién déi mét s6, mét bién hay mt d6i tugng nao do thanh kiéu 
Symbolic ta cé thé str dung m6t trong cac cach sau: 

>>s=sym(A) 

>>x=sym(x) 

>>syms x y 2% khai bao két hop x, y va z la bién symbolic 


Tinh dao ham bang ham diff cUla symbolic: Néu S 1a biéu thUtc symbolic 
thi: 


diff(S)dao ham cua S theo bién ty do 
diff(S,’v’)dao ham cla S theo bién v 
diff(S,’v’,n)dao ham cap n cla S theo v. 

Vi du: Tinh dao ham cua y = sinx3. 

>> syms x% khai bao x la bien kieu symbolic 
>> y=sin(x3); 

>> z=diff(y)% dao ham cua y 

z= 

3*cos(xA3)*xA2% sinh vien kiem tra ket qua 
>>pretty(z)% hien thi dang quen thuoc 

3 cos(x3) x2 


>>ezplot(x,y)% ve y theo x 


[missing_resource: graphics1.wmf] 
Hinh 3.1 — Vé d6 thi ham symbolic 
Tinh vi phan bang ham int - Néu S 1a biéu thtfc Symbolic thi: 


int(S)tich phan kh6ng xéc dinh cUa S theo bién mac nhién (muOn biét 
bién mac nhién nay ta ding ham findsym). 


int(S,v)tich phan khéng xac dinh cua S theo v. 
int(S,a,b)tich phan xac dinh cUa S trén can [a,b] 


int(S,v,a,b)tich phan xac dinh cua S theo v trén can [a,b] 


Vi du: Tinh j 2x 194122" gx 
>>syms X 
>>S=2*KN2*(19+12*xA2)/(7*(XA2+1)) 
>>y=int(S,x,0,1)% tich phan S theo x trén can [0,1] 
>>subs(y)% d6i sang kiéu s6 
Giai hé phuong trinh bang ham solve: 
>>help solve 
>>syms X Y 
>>[x,y]= solve('*x/2*sin(x42)-3*y=7','x+y=1') 


Sau khi thu duc nghiém x va y, sinh vién hay thay vao 2 phuong trinh 
trén va nhan xét két qua. 


Vé mat 3D bang ham ezsurf: Vi du vé mat S = f(z,y) = trén 


mién xac dinh: -5<x<5; -2 <y<2. 


¥ 
1+a22+y? 


>>syms X y 
>>S=y/(1+xA2+y2) 


>>ezsurf(S, [-5 5 —2*pi 2*pi]) 


Hinh 3.2 — Vé d6 thi 3D cho ham s6 symbolic 


Simulink 


Simulink 5.0 (Simulation and Link - R13) dugc MatWorks gidi thiéu vao 
thang 6 nadm 2003. No cho phép phan tich, m6 hinh héa va m6 phOng cac 
hé th6ng dong tuyén tinh va phi tuyén, lién tuc va rOi rac mOt cach truc 
quan trong m6i truOng giao tiép d6 hoa, bang cac thao tac chudt don 
gian. Co thé ndi, khéng tan dung du@c Simulink 14 mOt thiét thoi 16n cho 
nguOi lam cng tac m6 phOng! 


Kh6i d6ng Simulink bang m6t trong cdc cach sau: 


nhap: >>simulink 
hoac nhap chu6t vao trén menubar cUa Matlab 
Thu vién simulink hién ra nhu hinh 3.3: 


TruGc tién, sinh vién hay nhap chu6t vao cdc thanh cuOn cUa thu vién dé 
c6 cai nhin than thién vé simulink. 


TU day, dé c6 thé tao mé hinh bang simulink, hay: 


e nhap chudt vao biéu tuong cUa thu vién simulink 
e chon: File — New — Model trong Menu cUa thu vién Simulink 
e chon: File — New — Model trong cUfa sO 16nh cUa Matlab 


[“ Simulink Library Browser 
File Edit View Help 
C) ce ei Find | 
=) Bl Simulink 
«| Continuous 
~~ BH Discrete 
~~ Be] Functions & Tables Discrete 
oo Be] Math 
>| Nonlinear Functions & Tables 
2-| Signals & Systems 
Math 


>| Sources 
fe) GR CDMA Reference Blockset 


Klanline sr 


Hinh 3.3 — Ctra s6 chinh cua thu vién Simulink 


[=] untitled -/O) x! 


File Edit View Simulation Format Tools Help 


oom ff itt — eZ 


Hinh 3.4 — MOi truOng soan thao cUa Simulink 


Ctta sO nay (hinh 3.4) cho phép ta ‘nhap - kéo - tha’ vao tUing kh6i chtfc 
nang trong thu’ vién simulink. Vi du, dat vao day khOi ‘Sine Wave’ trong 
thu vién cUa (hinh 3.5): 


he 


15) x! 


Simulink Library Browser 
= Edit View Help 
B+} Continuou_& 


andon 
Num... 
Repeati 
Sequ... 


: >| Nonlinear Signal 
oo BH Signals & £ Gen... 


Ea ei CDMA Referen 
+) Bal Communication — 
‘yo Mak! Contenl Svete 

» 


[eluntitied *  aueTeT Es 
File Edit View Simulation 
Format Tools Help 


Sine Wave 


Hinh 3.5 — Lay m6t kh6i tl thu vién 

Sau khi dat tat ca cdc khOi can thiét cUla m6 hinh, ta nOi ching lai bang 
cach ‘nhap - git’ va kéo m6t duOng tU' ng6 ra cua kh6i nay dén ng6 vao 
cUa kh0i khac r6i tha phim trai chudt, m6t két ndi sé du@c thiét lap. 


XAy dung m6 hinh hé thong xe tai: 


by 
(friction) 


Hinh 3.6 — M6 hinh xe tai 
cho bdi phuOng trinh: m= = u — bv hay a = —(u—bv). 


Trong dé m 1a kh6i luOng xe, u 1a luc tac dOng cUa d6ng co (ngé vao cla 
m6 hinh), b 1a hé sO ma sat va v 1a van t6c dat dudc (ng ra cUa m6 hinh). 


Xuat phat diém cUa viéc xay dung cdc m6 hinh hé thOng tl cdc phuong 
trinh vi phan tuong Ung 1a cdc khOi tich phan (Integrator). Néu trong 
phuong trinh mé ta hé thOng cé vi phan bac n thi ta sé dat vao m6 hinh n 
khdi tich phan, do quan hé [ a“ =0, 


¢ MO m6t ctta s6 mé hinh mdi. 

e Dat vao m6 hinh khOi ‘Integrator’ tl’ thu vién ‘Continuous’ va ké cac 
duOng thang nOi dén ngé vao va ng6 ra cUa kh6i nay. 

e Dat nhan ‘vdot’ (dv/dt) cho cho dung ni dén ng6 vao va ‘v’ cho 
duOng nOi dén ng6 ra bang cach nhap dup chudt ngay phia trén cac 
duOng nay. 


vdot 4 v 
5 


Integrator 


TU phuong trinh hé théng ta thay dv/dt (vdot) bang tich cUa thanh phan 
(1/m) va thanh phan t6ng (u-bv), nén ta thém kh6i 1/m ngay tru6c kh6i 
tich phan: 

e Dat vao khOi ‘Gain’ trong thu vién 


>| Math 


e Nhdp dup chu6t vao khOi nay dé thay d6i d6 loi thanh 1/m. 
e Dat nhan ‘inertia’ cho khOi nay dé tugng trung cho quan tinh cUa xe 
(nhap dup vao nhan ‘Gain’ bén dui kh6i). 


wdot 4 v 
5 


inertia Integrator 


Bay gid ta dat kh6i tng vi 2 ngé vao ‘+-‘, ngd vao ‘+’ sé dudc ndi vi 
u, ngd vao ‘-’ sé dudc n6i vGi thanh phan bv dé dudc (u-bv) 


e Dat vao khOi ‘Sum’ trong thu vién 
2] Math 
e Nhdp dup vao kh6i nay dé d6i ngé vao tl ‘++’ sang ‘+-’ 
**EEKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Dé dudc thanh phan bv ta chi can dat thém kh6i ‘Gain’ vi d6 loi b: 


° Dat kh6i ‘Gain’ cd dé loib 
e Dat nhan 1a ‘damping’ tung trung cho thanh phan IUc can cUa xe. 


*EEKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Dén day viéc xay dung m6 hinh xe tai vGi ngé vao u va ng6 ra v coi nhu 
hoan thanh. Tuy nhién, dé m6 phong m6 hinh nay, ta can dat thém kh6i 


‘Step’ vao u va hién thi v trén khi ‘Scope’ 


t kh6i ‘Step’ trong thu vién bién dO u ngay ngé vao. 
t kh6i ‘Scope’ trong thu’ vién 


. Da 
. Da 

2 Sinks 
ngay ngo ra v. 

**EKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 

Nh6 rang m, b va u la cac bién can dugc gan tri truGc khi mé phong. 

>>m=1000 

>>b=50 

>>u=500 


ThOi gian m6 phong hé thOng tty thudc vao thong s6 ‘Stop time’ trong 
menu ‘Simulation Simulation paramrters’, gia sU' dat 120 (Hinh 3.7). 


Dé chay m6 phOng, ta co thé thuc hién bang 1 trong cdc cach: 
e nhap chu6t vao biéu tuong trén menubar cUa m6 hinh 
e chon: Simulation Start 
¢ Ctrl-T 


Nhap dup vao ‘Scope’ dé xem két qua m6 phong. 


Simulation Parameters: untitled 


Variable-step a4 ode46 (Dormand-Prince) 7” 


Lai | a= 
fave | Jao | 
[ave 


i 
ox | cancer| etn | avn | 


Hinh 3.7 — Thay d6i thong s6 m6 phong 


Xay dung m6 hinh hé th6ng diéu khién vi tri motor DC cho bdi phuong 
trinh vi phan sau: 


do. 1 . dd 
ae ee oe 


d_1 pg; _ do 


Trong do: 


J = 0.01 Kgm2/s2 14 moment quan tin cUa rotor 
b = 0.1 Mms 1a hé s6 ma sat cUa cac bO phan co khi 
K = Ke = Kt = 0.01 Nm/A la hang s6 stic dién dng 
R = 10 ohm 1a dién tré day quan 
L=0.5H1ahé sé tucam 
V 1a dién dp dat 1én cuOn day cua motor 

1a vi tri tr'uc quay (ng6 ra cUa m6 hinh) 


ila dong dién chay trong cuOn day cUa motor. 


ba 


Hinh 3.8 — M6 hinh toan hé diéu khién vi tri motor DC 


Quan sat tung phUOng trinh m6 ta hé thOng ta thay cau tric cUa ching 
citing tuong tl nhu phuong trinh cUa cau . Sinh vién lan luot thc hién 
trén tung phuOng trinh dé dugc (hinh 3.9): 


File Edit View Simulation Format Tools Help : 


| fc difdt B i 
+ _ 
= $ 


Sum Inductance Integrator 


i d2(thetaydt2 Ei d(thetaydt 


Sum1 Inertia Integrator1 


Fe 


(Hinh 3.9) 


Ké@t hop 2 phuong trinh: 


i=} motor (5) x} 


File Edit View Simulation Format Tools Help 


Sum Inductance Integrator 


<s 


damping 


Re [100% jode4s 7 


(Hinh 3.10) 


Dat vao m6 hinh kh6i ‘Step’ dé lam tin hiéu tham khao, khOi ‘Scope’ dé 
quan sat dap Ung. Sinh vién hay gan tri cho tat ca cdc thong s6 cua m6 
hinh, thu'c hién m6 phOng va quan sat dap Ung (Hinh 3.11). 


Hay thay d6i ng6 vao m6 hinh bang kh6i tao xung vung. M6 phong, 
quan sat két qua va nhan xét. 


TU chOn 


Sinh vién hay tinh dao ham cap 2 cla ham y = xe — bang tay va kiém 
chung két qua bang symbolic. 


CO 
Tinh tich phan sau va kiém chttng két qua bang symbolic: [ e~*” 
0 


[= motor (5) x! 


File Edit View Simulation Format Tools Help 


Resistance 


“> difdt i 


Inductance Integrator 


Coe eae 
Kt — = 
Sum1 Inertia Integrator1 


damping 


oom | fet — te 


Integrator2 Scope 


Hinh 3.11 - M6 hinh Simulink hoan chinh cla hé diéu khién vi tri motor 
DC 


Xay dung m6 hinh hé thOng xe IUta cho bdi phuong trinh: 


= F—k(a — 22) — pM, gS 


dz 
M ! 
1 dt2 


= k(21 — 22) — pMyg 


d2x 
M = 
2 dt 


r= a 


(Hinh 3.12 Photo courtesy: Dr. Howard Blackburn) 

Trong do cdc thong s6 tuOng trung nhu sau: 

M1=1 kg 1a khdi luOng toa kéo; 

M2=0.5 kg 1a khdi lugng toa khach; 

k=1 N/sec 1a d6 cUing 16 xo két n6i gilfa 2 toa; 

F=1 N 1a luc tac d6ng cUla dau may (ng6 vao m6 hinh); 
=0.002 sec/m 1a hé sO ma sat lan; 

g = 9.8 m/s/2 la gia tOc trong trudng 


x1, x2 vi tri 2 toa (ng6 ra). 


=. 
“Mgt, HM 8X5 


Hinh 3.13 — M6 hinh todn cUa hé thong xe IUta 


M6 hinh hé théng vién thong 
Day 1a gido trinh vé m6 hinh hé th6ng vién th6éng 
MUc tiéu 


Bai thi nghiém nay tim hiéu phuong phap m6 hinh héa m6t hé thong 
théng tin don gian nhu hinh 4.1: 


*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 

Hinh 4.1 - M6t m6 hinh truyén tin tiéu biéu 

Ching ta sé xt ly tin hiéu ngu6n (ma héa, kiém soat 16i, di€u ché) va 
truyén qua m6t kénh truyén cé nhiéu sau d6 khéi phuc tin hiéu nhan 
duQc (giai diéu ché, giai ma) va xac dinh x4c suat 106i tin hiéu. 
Tham khAo 


[1].The Mathworks Inc., Matlab User’s Guide — Communications Toolbox, 
2003. 


(2].Nguyén Httu Tinh - Lé Tan Hing - Pham Thi Ngoc Yén - Nguyén 
Thi Lan Huong, Co sO Matlab & Ung dung, NXB KH va Ky thuat, 1999. 


[3].http://murray.newcastle.edu.au/uers/staff/eemf/ELEC352/notes.htm 


ThUc hanh 


Dé cé thé thuc tap tt bai thi nghiém, sinh vién can cé kién thc vé XU 
ly sO tin hiéu (Digital Signal Processing), Truyén dif li€u (Data 
transmittion) va CO sO Vién thong. Do dé, bai thi nghiém nay khéng bat 
buGc d6i vdi sinh vién chuyén nganh Tin hOc (néu c6) va sinh vién Dién 
tt’ hudng DKTD. Trong trudng hop do, sinh vién cé thé chon bai 5 hodc 
bai 7 dé thu tap. 


Xay dung mo hinh hé thong thong tin ndi chung 1a mOt viéc rat 16n, vGi 
nhiéu van dé can gidi quyét. Do dé trong pham | vi bai thi nghiém nay 
chiing ta chi cé thé khao sat mOt s6 khau trong tOng thé cua mOt hé 
thong truyén thong. 


Hé th6ng thong tin lién tc (Analog Communications) 


Trong phan nay ching ta sé tim hiéu cdc phuong phap diéu ché va giai 
diéu ché tin hiéu lién tuc (Analog Modulation and Demodulation). 


Diéu ché bién dé séng mang bi nén hai bang canh (Double-Sideband 
Suppressed Carrier Amplitude Modulation — DSB-SC): Phuong phap diéu 
ché nay dUa theo céng thtic sau: 

y(t) = m(t)cos(2 fct+ c) 


Trong d6: m(t) 1a tin hiéu ha tan (théng tin) ding dé diéu ché séng mang 
tan sO cao c(t) = cos(2 fct+ c) dé thu dugc tin hiéu diéu ché y(t). 


*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 

Hinh 4.2 — Nguyén tac diéu ché bién d6 

Ta co thé thuc hién phuong phap diéu ché nay trong Matlab nhu sau 
(sinh vién nén ghi trong file .m dé thuan tién cho viéc su dung lai 6 cdc 
cau sau): 

>>Fc=10;% tan sO séng mang [Hz] 

>>Fs=60;% tan sO lady mau dé mé phong 

>>Fm=1; % tan s6 tin hiéu ha tan 

>>t=[0:1/Fs:3];% time window 


>>c=cos(2*pi*Fc*t);% tin hiéu s6ng mang 


>>m=sin(2*pi*Fm*t);% tin hiéu ha tan 1Hz 
>>subplot(311), plot(t,m)% vé tin hiéu ha tan 
>>subplot(312), plot(t,c)% vé song mang 

>>y=m.*c;% diéu ché 

>>subplot(312), plot(t,y)% vé tin hiéu diéu ché 
>>z=abs(fft(y));% ph6 tin hiéu diéu ché 
>>z=z(1:length(z)/2+1);% lay 1/2 
>>frq=[0:length(z)-1]*Fs/length(z)/2;% frequency window 
>>figure% new figure 

>>plot(frq,z);% vé pho 


Sinh vién hay quan sat ph6 tin hiéu diéu ché, va hay so sanh vGi tan s6 
song mang. 


Thu ding ham amod cua Communications toolbox thUc hién cdng viéc 
trén. 


>>help amod 
>>Y=amod(m, Fc, Fs, 'amdsb-sc’);% diéu ché 


>>subplot(313),plot(t, Y) 


[missing_resource: graphics1.wmf] 


Hinh 4.3 — Tin hiéu diéu ché bién d6 


TuOng tu nhu trén, sinh vién hay thuc hién diéu ché bién d6 séng mang 
bi nén don bang canh (Single-Sideband Suppressed Carrier Amplitude 
Modu-lation — SSB-SC). So sdnh v6i phuOong phap diéu ché bién d6 DSB- 
no OF 

Giai diéu ché bién d6: C6 nhiéu phuong phap gidi diéu ché bién dd. O 
day ta xét phuOng phap giai diéu ché d6ng b6 (Synchronous 
Demodulation), 4p dung cho phuong phéap diéu ché bién d6 vGi séng 
mang bi nén SC (Suppressed Carrier) va phuong phap diéu ché bién dd 
vGi song mang duc truyén TC (Transmission Carrier), theo m6 hinh sau: 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 

Hinh 4.4 — Nguyén tac giai diéu ché bién dé 


Thuc hién giai diéu ché nhu sau (phai dam bao rang tin hiéu diéu ché y 
con trong WorkSpace cUa Matlab): 


>>y_mix=y.*c; % trOn tin hiéu diéu ché vi séng mang 
>>[num,den]=butter(3, (Fc/2)/(Fs/2)); % low-pass filter tan s6 cat Fc/2 
>>m_r=filter(num,den,y_mix); % cho tin hiéu da trOn qua b6 loc 
>>close all% 
>>plot(t,sm)% vé tin hiéu ha tan ban dau 
>>hold on% 
>>plot(t,m_r,'r')% vé tin hiéu ha tan vlra hi phuc 

Gidi thich vi sao ta ding bO lOc thong thap vGi tan s6 cat Fc/2. 


Sinh vién hay quan sat va so sanh su khac biét gitfa tin hiéu giai diéu 
ché va tin hiéu ha tan ban dau. Giai thich. 


Thtr ding ham ademod cUa Communications toolbox dé gidi diéu ché: 


>>help ademod 
>>M_R=ademod(y,Fc,Fs,'amdsb-sc');% giai diéu ché 


>>plot(t,M_R,'g’) 


Hé th6ng thong tin rdi rac (Digital Communications) 


Trong phan nay, ching ta sé thc hién m6 hinh m6t hé thOng truyén 
thong rdi rac nhu so dé sau: 


**EEKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 4.5 — M6 hinh m6t hé th6ng truyén tin rdi rac 
Sinh vién can doc phan help ca tng 1@nh hiéu va biét cach st? dung. 
Tao tin hiéu nguOn: 
>>clear 
>>close all 
>>K=11;% chiéu dai tU tin (message length) 
>>msg=randint(K*2,1);% tao tin hiéu nguOn (digital message) 
>>subplot(411),stem(msg,’.’)% vé tin hiéu ngudn 
>>ylabel(‘digital message’) 
Ma hoa BCH: 
>>N=15;%chiéu dai tu’ ma (codeword length) 


>>code = encode(msg,N,K,'bch');%ma héa BCH 


>>subplot(412), stem(code,’.’) 
>>ylabel(‘BCH’) 


Diéu ché QASK: SU dung b6 diéu ché s6 M-QASK vGi M=16 (sinh 
vién can xem lai ly thuyét vé bo diéu ché nay) 


>>M=16;% multiple number 

>>Fd=1;% t6c d6 nguén tin 

>>Fc=10;% tan s6 séng mang 

>>Fs=30;% tan s6 lady mau 

>>modu = dmod(code,Fc,Fd,Fs,'gask',M);% di€u ch6 M-QASK 
>>subplot(413),plot(modu)% vé tin hiéu diéu ché 
>>ylabel(‘QASK’) 


Kénh truyén co nhiéu: Tin hi€éu sau khi diéu ché duc truyén trén kénh 
truyén co nhiéu trang Gaussian c6ng suat 0.1W. 


>>std_value=0.1;% cong suat nhieu 

>>modu_noise = modu+randn(length(modu),1) std_value; % add noise 
Giai diéu ché QASK: 

>>demo = ddemod(modu_noise,Fc,Fd,Fs,'gask', M);% giai diéu ché 
Giai ma BCH: 

>>msg_r = decode(demo,N,K,'bch');% tin hiéu thu 

>>subplot(414), stem(msg_r,’.’)% vé tin hiéu s6 thu dudc 


>>ylabel(‘received message’) 


[missing_resource: graphics2.wmf] 
Hinh 4.6 — So saénh tin hiéu truyén/nhan 

Xac dinh 16i: 
>>rate=biterr(msg,msg_r,M) 

Sinh vién hay ding ham channel cUa Communications toolbox dé m6 hinh 
héa kénh truyén 6 buGc_ . Lap lai cdc qua trinh giai diéu ché va giai ma, 
so sanh két qua (lUu y: bang thong cUa kénh truyén phai du 16n). 

TU chOn 


Thuc hién diéu ché va giai diéu ché tan s6 FM vGi Fm, Fc, Fs 6 cau 
III.1. . 


Thuc hién m6 hinh Digital Communications nhu hinh sau: 


Gidi ma 


msg tf 
1 Ham 


ming 


M6 hinh hé théng diéu khién tu déng 
Day la gido trinh vé m6 hinh hé thOng diéu khién tu déng 
MUc tiéu 


Trong bai thi nghiém nay ching ta sé tim hiéu phuong phap m6 hinh hdéa m6t hé 
diéu khién tu dOng, bao g6m: 


e Ham truyén va phuoOng trinh trang thai cUa hé thong 

e Dap Ung vong ho va dap Ung vong kin cua hé thong 

e Xay dung b6 diéu khién PID 

¢ Chinh dinh théng s6 cUa b6 diéu khién va khao sat dap Ung cUa hé thOng. 
***ESORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


Hinh 5.1 — M6t m6 hinh hé théng diéu khién tiéu biéu 


Tham khAo 
[1].The Mathworks Inc., Matlab Notebook User’s Guide — Control toolbox, 2003. 


[2]. Pham Van Tan, Bai gidng m6n Co sé Tu dong hoc, BO mén Vién Thong va 
Tu d6ng héa, khoa Cong nghé Thong tin, Dai hoc Can Tho, 2001. 


[3].Nguyén Céng Dinh, Phan tich va T6ng hp cac hé théng Diéu khién bang 
may tinh, NXB Khoa hoc va KY thuat, 2002. 


[4].http://www.engin.umich.edu/group/ctm 


[5]. http://www.shu.ac.uk/schools/eng/teaching/rw/pidtutorial.htm 


ThUc hanh 


Dé c6 thé thurc hién tt bai thi nghiém, sinh vién can nam ving cac kién thc 
co ban vé Diéu khién tu d6ng (Co sé tu déng hdc). Do do, bai nay khéng bat 
buGc d6i vGi cac sinh vién Tin hc (néu c6) va cac sinh vién Dién tl theo hudng 
Vién théng. Trong truOng hgp dé, co thé sinh vién da thUc tap bai 4 hoac sinh 
vién c6 thé chuyén sang bai 7. 


Ham truyén va phuOng trinh trang thai cUa hé thong 

Trong diéu khién tu’ dong, ngudi ta thuOng biéu dién mét hé th6ng vat ly bang 
ham truyén (transfer function) hay phuOng trinh trang thai (state-space equation) 
cUa n6 (di vGi cdc hé phi tuyén, dé dat duoc diéu nay, ngudi ta phai ding 
phuong phap tuyén tinh hdéa tung doan). 

Gia sU c6 hé th6ng diéu khién t6c d6 motor DC nhu hinh vé 5.2 [4]. Trong do: 
J = 0.01 kgm2/s2 14 moment quan tinh cUa rotor 
b= 0.1 Nms hé s6 ma sat 
K=Ke=Kt=0.01 Nm/Amp cac hang s6 sttc dién d6ng 
R= 1 ohm dién trd 
L=0.5H dién cam 
I:dong dién chay trong cuOn day cUa motor 
V: dién ap trén hai dau cuOn day motor — ngé vao 

: vi tri truc — ng6 ra 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 5.2 — M6 hinh todn m6t hé diéu khién t6c d6 motor DC 
Phuong trinh vi phan m6 ta hé th6ng nhu sau: 

J +04 =Ki 

L&+Ri=V-Ke 

Ham truyén: Bién d6i Laplace 2 vé cUa phuOng trinh trén ta dugc: 


s(Js + b)O(s) = KI(s) 
(Ls + R)I(s) = V — KsO@(s) 


Suy ra: [(Ls + R)(Js + 6) + K?]s0 = KV hay % = qpGermr? 


Biéu dién ham truyén nay trong Matlab ta thc hién nhu sau (sinh vién nén luu 
thanh file.m): 


>>J=0.01; 
>>b=0.1; 
>>K=0.01; 
>>R=1; 
>>L=0.5; 
>>num=K;% tU’'s6 cUa ham truyén 
>>den=[(J*L) ((J*R)+(L*b)) ((b*R)+K/2)];% mau sO ham truyén 
>>hamtruyen = tf(num,den) 
Dap Ung buGc vong hd: 
>>step(num,den)% hoac 


>>step(hamtruyen) 


[missing_resource: graphics1.wmf] 
Dap Ung xung vong he: 
***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


>>impulse(hamtruyen) 


[missing_resource: graphics2.wmf] 
Phuong trinh trang thai: Dang tOng quat: 


X =AX+BU 
Y= CX DU 


vGi X 1a vécto trang thai, U 1a vécto tin hiéu vao va Y 1a vécto tin hiéu ra. 


Bién trang thai va phuOng trinh trang thai: TU phuOng trinh vi phan m6 ta hé 
thong, néu dat x; = Ova x2 =i, ta CO: 


eee aac ae = eae x+| |v 
L2= ~4£ x, - ary + LV hay — —# + 
Yjar- 072 ov Y=(1 0X 


Biéu dién phuong trinh trang thai trong Matlab nhu sau: 
>>J=0.01; 
>>b=0.1; 
>>K=0.01; 
>>R=1; 
>>L=0.5; 
>>A = [-b/J K/J; -K/L -R/L]; 
>>B = [0;1/L]; 
>>C=[10]; 
>>D=0; 

Dap Ung bu‘6c vong ho: 
>>step(A,B,C,D) 

Dap Ung xung vong hd: 
>>impulse(A,B,C,D) 


Ta co thé chuyén d6i qua lai giffa ham truyén va phuoOng trinh trang thai bang 
16nh sau: 


>>[num,den]=ss2tf(A,B,C,D)% tlr PT trang thai sang ham truyén 


>>[A,B,C,D]=tf2ss(num,den)% tlr ham truyén sang PT trang thai 
Khao sat dap Ung vong hd cUa hé thOng d6i vi tin hiéu bat ky 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 5.3) 
Phai dam bao rang trong Workspace con bién hamtruyen cUa cau_, sinh vién c6é 
thé dung 1énh Isim dé khao sat dap Ung cUa hé d6i v6i tin hiéu bat ky. Gia sUr 
do 1a tin hiéu sin: 
>>close all 
>>t=0:0.1:2*pi; 
>>u=sin(pi/4*t); 


>>I|sim(hamtruyen,u,t)% mo phong dap ung voi tin hieu vao u 


[missing_resource: graphics3.wmf] 


BO diéu khién PID 

Cau tric m6t hé th6ng diéu khién PID nhu hinh sau: 
***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 5.4 — So d6 kh6i hé diéu khién PID 


Kps*+Kps+Ky, 
8 


Trong d6 ham truyén cUa khau PID la: Kp + Ar + Kps= 
v6i: KP 1a d6 loi cUa khau ti 16 (Proportional gain) 
KI 1a d6 oi cUa khau tich phan (Integral gain) 


KD 1a d6 Igi khau vi phan (Derivative gain) 


Viéc hiéu chinh phi hgp 3 théng s6 KP, KI va KD sé lam tang chat lugng diéu 
khién. Anh hung cla 3 thong sO nay lén hé th6ng nhu sau: 


ae T. gian 

ae. (||| Manse 
time) 

response) 

KP Giam 

KI Giam 

man It thay 
déi 


Vot 
16(Overshoot) 


Ta ng 


Tang 


[missing_resource: graphics4.wmf] 


BO diéu khién ti 16 P: 


T. gian qua 
d6(Settling 
time) 


It thay d61 


Tang 


Giam 


EE SORRY, THIS MEDIA TYPE IS NOT SUPPORTED,*** 


Hinh 5.5 — BO diéu khién ti 16 P 


Sai sO xac 
lap(Steady- 
State err.) 


Ta ng 


Khong xac 
dinh 


Thay d6i it 


Thuc hién trong Matlab: Ta cé ham truyén cUa motor DC nhuIII.1. : 


>>J=0.01; 
>>b=0.1; 
>>K=0.01; 


>>R=1; 


>>L=0.5; 

>>num=K; 

>>den=[(J*L) ((J*R)+(L*b)) ((b*R)+K‘2)]; 

Khi thém vao khau ti 16 P, ta cd ham truyén vong hé: 

*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
>>Kp=100; 

>>numa=Kp*num; 

>>dena=den; 

Xac dinh ham truyén vong kin cUa hé thOng ta ding 1énh cloop: 
>>[numac,denac]=cloop(numa,dena) 

Dap ting Step vong kin ca b6 diéu khién ti 16 nhu sau: 
>>t=0:0.01:2; 


>>step(numac,denac) 


[missing_resource: graphics5.wmf] 


Sinh vién hay so sénh vGi dap Ung cUa hé khi chua c6 b6 diéu khién ti 16, 6 
cau III.1. (luu y dén cac thing s6: thdi gian lén, d6 vot 16, thi gian qua dd). 


TuOng tU, sinh vién hay so sénh vGi dap Ung xung. 

BO diéu khién Vi tich phan ti 16 PID: 
***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 5.6 — BO diéu khién PID 


Khi thém b6 diéu khién PID, ham truyén h@ cla hé th6éng 1a: 


*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 

>>Kp=100; 

>>Ki=1; 

>>Kd=1; 

>>numc=[Kd, Kp, Ki]; 

>>denc=[1 0]; 

>>numa=conv(num,numc);% tich chap ttr s6 

>>dena=conv(den,denc);% tich chap mau sO 

Ham truyén vong kin hi ti€p 4m don vi: 

>>[numac,denac]=cloop(numa,dena); 

Dap Ung Step cUa hé diéu khién PID: 

>>step(numac,denac) 
Sinh vién so sénh vGi dap Ung cUa b6 diéu khién ti 16 P G cau _, nhan xét. 
Dua vao bang tOng két anh huGng cUa KP, KD va KI d6i vGi hé thong diéu 


khién, sinh vién hay thay d6i 3 thong s6 nay va kiém chting dap Ung cla hé 
thong. 


[missing_resource: graphics6.wmf] 


Hiéu chinh thong s6 cla b6 diéu khién PID 

M6t phuong phap c6 dién nhung don gian va hiéu qua dé chinh dinh 3 théng sO 
KP, KI va KD cla b6 diéu khién PID 1a phuong phap Ziegler-Nichols (Ziegler 
Nichols Tuning Method). Thu tuc chinh dinh nhu sau: 


1. Chi diéu khién hé th6ng bang bO diéu khién ti 16 KP (dat KI=KD=0). 


2. Tang KP dén gia tri KC ma 6 dé hé thong bat dau bat 6n (bat dau xuat hién 
su giao dOng - diém cu cUa ham truyén kin nam trén truc ao j ). Xac dinh tan 
sO ccUa giao dOng vtra dat. 


TU 2 gid tri KC va _c vUa dat, cdc thong s6 sO KP, KI va KD dugc xac dinh nhu 
bang sau: 


BO diéu khién KP KI KD 

P (ti 16) 0.5 KC 

PI (tich phan ti 14) 0.45 KC 0.191KP c 

PID (vi tich phan ti 16) 0.6 KC 0.318KP c 0.785KP/ c 
PID vGi mt it vot 16 0.33 KC 0.318KP c 2.07KP/ c 
PID khéng vot 16 0.2 KC 0.53KP c 3.14KP/ c 


3. Tinh chinh lai 3 thong s6 nay dé dat duc dap Ung nhu mong muOn. 
Vi du: Gia st? can thiét ké b6 diéu khién PID cho hé th6ng sau: 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Buéc 1: Diéu khién hé thOng chi vGi b6 diéu khién ti 1é: 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Buc 2: Xac dinh KC va _c ma 6 do hé thOng bat dau giao d6ng - ding ham 
rlocus cUa Matlab (sinh vién nén lutu thanh file .m hodc thao tac trong Matlab 


Editor sau d6 copy va dan vao Workspace ca doan lénh dé dé dang cho viéc 
hiéu chinh cdc thong s6 6 phan sau): 


>>close all 


>>num=5; 

>>den=[1 10 100 0]; 
>>[numc,denc]=cloop(num,den); 
>>htkin=tf(numc,denc)% ham truyen vong kin 
>>rlocus(htkin);%ve qui dao nghiem 


>> axis([-10 10 -15 15]) 


[missing_resource: graphics7.wmf] 
Xac dinh Kc va c bang ham rlocfind: 
>>[Kc,Omegac] = rlocfind(htkin) 


Nhap chu6t vao diém giao nhau gitta qui dao nghiém va trUc ao cUa dé thi, 
trong WorkSpace ta dugc: 


Kc = 

199.5793 

Omegac = 

-10.0145 

0.0072 +10.0072i 

0.0072 - 10.0072i 

Nhu vay ta dug@c KC=200 va c = 10. Suy ra thong s6 cUa b6 diéu khién PID: 
KP = 0.6KC = 120 

KI = 0.318KP c = 381.6 


KD = 0.785KP/ c = 9.4 


Thu dap Ung cUa hé: 

>>Kp=120; Ki=381.5; Kd=9.4; 

>>numc=[Kd, Kp, Ki]; 

>>denc=[1 0];% ham truyen cua PID 
>>[numac,denac]=cloop(conv(num,numc),conv(den,denc)) 
>>step(numac,denac) 


BuGc 3: Thu hién tuong tu nhu'III.2. , sinh vién hay diéu chinh m6t lugng 
nho 3 thong sO KP, KD va KI dé dudc dap Ung t6t hon. 


[missing_resource: graphics8.wmf] 
Sinh vién hay thiét ké b6 diéu khién PID cho hé thOng sau: 


***E SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


TU chOn 


Sinh vién hay thiét ké b6 diéu khién Vi phan ti 16 (Proportional-Derivative 
controller): 


***E SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


Sinh vién hay thiét ké b6 diéu khién Tich phan ti 16 (Proportional-Integral 
controller): 


***E SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


Tao giao dién trong Mathlab 
Day 1a gido trinh vé tao giao dién trong Matlab 


MUc tiéu 

Bai thi nghiém nay trang bi cho sinh vién mOt sO kién thttc co ban dé xay 
dung giao dién nguodi ding trong méi truOng Matlab, nham hoan thién 
mOt chuOng trinh Ung dung nhat dinh. 

Tham khAo 

[1].The Mathworks Inc., Matlab Notebook User’s Guide, 2003. 


[2].Nguyén HUtu Tinh - Lé Tan Hing - Pham Thi Ngoc Yén - Nguyén 
Thi Lan Huong, Co sO Matlab & Ung dung, NXB KH va ky thuat, 1999. 


[3].Nguyén Hoai SOn - DO Thanh Viét - Bui Xuan Lam, Ung dung 
MATLAB trong tinh todn ky thuat, Tap 1, NXB DHQG Tp. HCM, 2000 
(trang 328-344). 


ThUc hanh 


Ciing nhu cac ng6n ngff cap cao khac, Matlab h6 tro nhiéu céng cu chtfc 
nang cho phép lap trinh tao giao dién su dung dep va nhanh chong. Vi 
du, cdc dang ntit an, cUfa sO soan thao, cdc dang menu, ... nhu hinh 6.1. 


***SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 6.1 — Cac cong cu h6 tro giao dién 


Trong bai thi nghiém nay, sinh vién lan luot thuc hién cdc buGc dudc néu 
ra dé dudc m6t giao dién don gian nhu hinh vé 6.2. Trong mOi bu6éc 
thuc hanh, sinh vién hay quan sat ky dap Ung cUa chuOng trinh, tu dé tu 
minh rut ra kinh nghiém vé viéc diéu chinh cdc thu6éc tinh d6 hoa cua 
Matlab. 


Tao cUa sO chinh figure — ‘Welcome to User Interface’s Giude’ theo cach 
sau: 


e M6 Matlab Editor 
e Nhdap ni dung sau: 


% User Interface's Guide 
% 
% Matlab Experiments 2003 


% TcAD, CIT, Cantho University 


% Initialize whole figure... 

namefig = 'Welcome to User Interface"s Guide’; 

figpos = get(0,'DefaultFigurePosition'); % lay vi tri mac nhien 
figpos(1)= figpos(1)-10; figpos(2)= figpos(2)-10; 

figpos(3)= figpos(3)+10; figpos(4)= figpos(4)+10; 

% Tao figure 

fig=figure( ... 

‘Name’, namefig, ... 

'NumberTitle’,'off’, ... 


‘Position’ ,figpos); 


Welcome to User Interface’s Guide aD x) 


File Edit Yiew Insert Tools Window Help 


Matlab Experiments 


Hinh 6.2 — M6t giao dién don gian 


e Luu thanh file giaodien.m 
e >>giaodien < >% thi hanh 


Quan sat két qua trén hinh 6.3 (Luu y cac thu6c tinh: Name, Position) 
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(Hinh 6.3) 


° >>set(fig)% xem thudc tinh cUa figure trong cUra s6 1énh 


Thém vao figure mOt axes cho phép hién thi d6 hoa: 


e Nhap ti€p n6i dung sau vao file giaodien.m 


% main axes 
axs=axes('Position’,[0.05 0.4 0.65 0.55)); 


e Luu va thi hanh file nay. 
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(Hinh 6.4) 


Hién thi dong ‘Matlab Experiments’ bén duGi axes: 


e Nhap tiép n6i dung sau vao file giaodien.m 


% text 
txtpos=[10 50 425 50]; 
txt=uicontrol(... 


‘Style’, 'text’,... 


‘BackgroundColor',[0.8 0.8 0.8]... 
'ForegroundColor'’,[0.4 0.5 0.3]... 
'String’,'Matlab Experiments’... 
‘Position’ ,txtpos,... 
"Fontname'’,'Courier’,... 
'FontWeight','Bold'’,... 
'FontSize',26); 


e Luu va thi hanh file nay. 
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(Hinh 6.5) 
Tao m6t frame c6 shadow bén phai figure dé dat cdc ntit chic nang: 


e Nhap ti€p n6i dung sau vao file giaodien.m 


% Console frames 


p1=0.755; p2=0.05; p3=0.2; p4=0.90; 
frm1ipos = [p1 p2 p3 p4]; 

frm2pos = [p1-0.005 p2+0.005 p3 p4]; 
% shadow frame 

frm1=uicontrol( ... 

'Style','frame’, ... 

'Units', ‘normalized’, ... 
‘Position',frm1pos, ... 
'ForegroundColor',[0.4 0.4 0.4],... 
'‘BackgroundColor',[0.4 0.4 0.4]); 

% main frame 

frm2=uicontrol( ... 

‘Style’,'frame’, ... 

"Units', "normalized, ... 
‘Position',frm2pos, ... 
'ForegroundColor',[0.7 0.7 0.7],... 
'‘BackgroundColor’,[0.65 0.65 0.65]); 


e Luu va thi hanh file nay. 
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(Hinh 6.6) 
Tao nit ‘Close’ c6 chUc nang dong cla s6 figure hién hanh: 


e Nhap tiép ndi dung sau vao file giaodien.m 


% Close button 

closeHndl=uicontrol(... 

‘Style',‘pushbutton’, ... 

"Units',‘normalized, ... 

‘Position',[p1+0.01 p2+0.05 p3-0.025 0.05], ... 
‘String’,'Close’, ... 

'Foregroundcolor','b’,... 

'Fontsize',9,... 


'Callback','close'); 


e Luu va thi hanh file nay. 
Van dé quan trong nhat d6i vGi m6t nut chUc nang 1a thi hanh cong viéc 
tuOng Ung khi nguOi su dung thao tac. Thu6c tinh ‘CallBack’ cho phép: 
thi hanh mt lénh cUa Matlab dui dang chu0i (bao g6m 1énh goi ham, 
script file, biéu thUc toan, ...). Sinh vién thUr thay 1€nh close bang mét 
1énh khac, chang han demos va quan sat dap Ung. 
Tao nut Sphere dé vé m6t hinh cau 3D: 


e Nhap tiép ni dung sau vao file giaodien.m 


% Sphere button 

sph=uicontrol( ... 

'Style',‘pushbutton’, ... 

"Units', "normalized, ... 

‘Position',[p1+0.01 p4-0.05 p3-0.025 0.05], ... 
'String’,'Sphere’, ... 

'Foregroundcolor’,'b’,... 

'Fontsize',9,... 

‘Callback',['[x,y,z]=sphere(35);’, 'surfl(x,y,z);', 'axis equal']); 


e Luu va thi hanh file nay. 
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(Hinh 6.7) 
Tao nit Picture dé hién thi anh mau: 


e Nhap tiép ndi dung sau vao file giaodien.m 


% Picture button 

pic=uicontrol( ... 

‘Style',‘pushbutton’, ... 

"'Units',‘normalized, ... 

‘Position’,[p1+0.01 p4-0.125 p3-0.025 0.05], ... 
‘String’,'Picture’, ... 

'Foregroundcolor','b’,... 

'Fontsize',9,... 


‘Callback’, ['imshow(imread("flowers.tif"))']); 


e Lut va thi hanh file nay. 
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(Hinh 6.8) 


Sinh vién luu y cach biéu dién nhiéu lénh 6 dang chudi cho ‘CallBack’. 
Trong truOng hOp co qua nhiéu lénh phuc vu chUc nang nay, ta nén dua 
ching vao m6t script file hoAc mt ham khac. Ngoai ra, néu ta tao giao 
dién dui dang mdt ham (function) thi thu6c tinh ‘CallBack’ cho phép goi 
md0t ham con du@c viét ngay trong file nay. 


Sinh vién hay tao thém ntit ‘About’ sao cho khi ngu0i sf dung An nit sé 
m6 cla sO mGi nam ngay gitfa man hinh 6.9. 


WeiCome fo 


(Hinh 6.9) 


TU chOn 


Sinh vién hay xem m6t sO file tao giao dién demo cUa cac toolbox trong 
Matlab chang han: graf3d.m trong \toolbox\matlab\demos, dé tham khao 
cach xay dung giao dién theo dang ham. VGi cach nay, thu6c tinh 
‘CallBack’ trong ‘uicontrol’ cho phép goi trUc ti€p mOt ham ma ham dé 
duQc xay dung ngay trong file giao dién (nghia la khong can tao m6t file 
* m khac phuc vu cho ‘CallBack’). 


Thu sw lai giao dién cUa bai thi nghiém nay theo dang trén. 


Thiét ké - m6 phong mach va vé mach in 
Day 1a gido trinh vé thiét ké - m6 phéng mach va vé mach in 


MUc tiéu 


Gitip sinh vién buGc dau lam quen vGi céc phan mém sau: Multisim 
V6.20 (Electronics WorkBench) ding dé mé phéng cdc mach tuOng tu va 
mach s6; OrCAD Release 9 diing dé vé Schematic, m6 phOng va vé 
mach in (PCB). 


Tham khAo 
[1]. OrCAD® Inc., User’s Guide - Capture, 1998. 
[2]. OrCAD® Inc., User’s Guide - Layout, 1998. 


[3].http://www.electronicsworkbench.com/tutorial/ 


[4].http://www.rit.edu/~labeee/ 


ThUc hanh 


Mac du Multisim va OrCAD/Pspice la cdc phan mém rat 16n vGi nhiéu 
tinh nang 4p dung manh va phuc tap, nhung sinh vién chuyén nganh Dién 
tUr hoan toan cé thé tu nghién cUtu dudc. Do dé, bai thi nghiém nay 1a 
khéng bat buéc d6i vGi tat ca sinh vién (dac biét 1a sinh vién Tin hoc, 
néu co). 


Trong truOng hgp sinh vién khéng chon bai thi nghiém nay thi c6 thé 
chon Control System Toolbox (d6i vGi sinh vién DKTD) hoac 
Communications Toolbox (d6i vGi sinh vién Vién Thong) hodc Image 
Processing Toolbox (d6i vGi sinh vién Tin hoc) cua MATLAB dé khao 
sat. 


Bai thi nghi€m nay tung d6i dai, sinh vién can sap xép th0i gian hgp ly 
dé thuc hién cac cong viéc sau: 


Multisim 


Trong phan nay, ching ta sé lam quen vGi Multisim 6.20 (hién da cé 
Multisim 2001 hoan thién hon), day 14 Version mi cUa Electronics 
WorkBench 5.x, né cho phép str dung két hop thu vién linh kién Analog 
va Digital vGi nhiéu cong cu than thién, h6 trd tuyét voi cho viéc m6 
phong cdc mach dién tu. 


Khoi d6ng Multisim: 


e Chay file Multisim.exe tl’ Start menu cua Windows 
e Nhdap vao biéu tuOng trén Desktop. 


Sinh vién hay nhap chu6t vao cac thu vién linh kién cUa Multisim trén 
Menu bén trai cla cUta sO chinh dé cé cai nhin t6ng quat vé cdc linh kién 
chUc nang. 


e Khi mu6n dat linh ki€n nao vao kh6ng gian thiét ké, chi can nhap 
chudt vao linh kién dé, r6i chon thong s6 phi hgp va dat vao vi tri 
mong muOn. 

¢ Dé tao mOt két n6i gitta 2 diém, chi can nhdp chudt vao diém tht? 
nhat r6i kéo dén diém tht? hai va nhap chu6ét lan nUia. 

*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 7.1 — Giao dién chinh cua MultiSIM 
Sinh vién hay vé va m6 phong mach dién sau: 
*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
Hinh 7.2 — MOt mach dién dién hinh 
e Ding Osilsoscope quan sat dang tin hiéu 6 chan 2-6 va chan 3 cua 


LM555. 
e Hay uGc lung tan sO giao d6ng va chu ky thao tac cUa mach. 


Kinh nghiém: D6i lic Multisim dua ra bang thong bao 106i trong qua trinh 
mo phong hay hiéu chinh th6ng s6 cUa cac linh kién. Truong hop do, ta 
g0 bo khoi Osilloscope, hiéu chinh va m6 phOng vai lan ri gan 
Osilloscope tré lai dé quan sat dang tin hi€u. Ngoai ra, néu trong mach 
c6 slr dung cac IC s6 (74xxx, 40xxx, ...) ta phai cé nguén Vcec=5V nhu 
hinh vé trén. 


OrCAD 


Cac chic nang chinh cua OrCad g6m: Vé so d6 mach (OrCAD Capture), 
m6 phong mach (OrCAD Pspice) va vé mach in (OrCAD Layout). Tuy 
nhién, trong phan thi nghiém nay ching ta khong dé cap dén OrCAD 
Pspice (day 1a mt trong nhttng chic nang rat manh cla OrCAD phién 
ban 9.0 trd lén, dé nghi sinh vién ty nghién ctu). 


Vé so d6 mach (Schematic) 
Kh6di dOng OrCAD Capture bang m6t trong cac cach sau: 


e Chon Capture CIS tur Start Menu cua Window 
e Nhdap vao biéu tuOng trén desktop 


Dé tao mOt Project mdi cho viéc vé so d6 mach ta thUc hién: 


e Chon File New Project... hodc nhap vao icon trén Menu bar 

e Khi hOp thoai ‘New Project’ xuat hién: Chon ‘Schematic’ va dat tén 
project trong menu ‘Name’ (gia st’ 1a BaiTN7). Co thé chon thu muc 
lu project nay trong menu ‘Location’. 


MOi truOng vé mach xuat hién khi ta nhdp dtip vao ‘PAGE1’ (hinh 7.3). 
TU day ta lan luot dat cdc linh kién vao va két noi dé dudc so dO mong 
muOn. 


i$ OrCAD Capture 


H-G Design Cache 
Lane Library 

so {) Outputs 

7 Referenced Prnierte 


| lo items selected Scale=100% X=30.48 Y=40.64 


Hinh 7.3 — MOi truOng vé mach dién cua OrcAD 

SU dung Menu tat sau: 

**EKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 7.4) 


Sinh vién hay vé so d6 mach nhu hinh sau. Trong dé cac thu vién linh 
kién dudc sl dung nhu sau: 
e Discrete: R, Capacitor Pol, Capacitor Non-Pol, Resistor Var, Photo 
PNP va Speaker. 
e Opamp: TL082 
e Amplifier: LM386 


Luu Schematic vUa vé vGi tén BaiTN7. 


[missing_resource: graphics2.wmf] 


(Hinh 7.5) 


Tao Netlist 


Sau khi luu so d6 vUa vé, dé tao Netlist phUc vu cho viéc vé mach in ta 
thu’c hién nhu sau: 


e Dong cUfa sO ‘Schematicl:Pagel’, quay vé Project Management 
‘BaiTN7’ 


etESORRY, THIS MEDIA TYPE IS NOT SUPPORTED :*** 
(Hinh 7.6) 


¢ Chon Tools — Annotate... dé danh s6 tht? ty’ cho linh kién. 

¢ Chon Tool — Create Netlist.... dé tao Netlist. 

e Khi cUfa s6 ‘Create Netlist’ xudt hién, ta chon layout Tab (file d& liéu 
netlist la BaiTN7.MNL), OK. 


Néu dam bao rang cé file BaiTN7.MNL trong thu muc hién hanh thi 
thoat khdi OrCAD Capture dé chudn bi cho viéc vé mach in. 


Luu y: Trong phan trén ching ta khong dé cap dén chic nang Kiém tra 
mach cUa OrCAD va m6t s6 chUtc nang h6 tro khac. 


Vé mach in (Printed Circuit Board-PCB) 
Khdi dOng OrCAD Layout bang m6t trong cac cach sau: 


e Chon Layout Plus trong Start Menu cua Window 
e Nhdap vao biéu tuOng trén desktop. 
e Chon File New dé tao m6t PCB layout m6i 


Sau thao tac nay, cla s6 ‘Load Template File’ xuat hign, ta can dat dung 
dan chi dén file DEFAULT.tch trong OrCAD/Layout_plus/DATA. Sau dé 
thc hién nap file Netlist vUia tao 6 buGc III.3.2 (vi du: BaiTN7.MNL). 


Dat footprint cho cac linh kién: 


Néu trong qua trinh vé Schematic ta khong chon footprint cho cac linh 
kién mOi khi dat vao thi 6 buGc nay OrCAD sé yéu cau chon kiéu chan 
cho cdc linh kién theo menu sau: 


Link Footprint to Component Wx] 


AutoECO cannot find a footprint for component Rl 
from part name R. 


Please choose one of the options below: 


(Hinh 7.7) 
Nhap vao menu ‘Link existing footprint to component ...’, chon Ia kiéu 
chan phi hQp cho tUing loai linh kién. O buc nay, kiéu chan linh kién 


can chon ding kich thuGc (c6 thé phai tao cdc footprint mdi). DOi vGi 
bai thi nghiém nay, sinh vién cé thé tham khdao bang sau: 


Linh kién Footprint Thu vién 
R Jumper600 JUMPER 


C CPCYL1/D.200/LS.100/.031 TM_CAP_P 


Photo NPN TO92 TO 
TL082, LM386 DIP.100/8/W.300/LS.100/.425 DIP100T 
VR VRES1 VRES 


Speaker CYL/D.150/LS.100/.31 TM_CYLND 


Sap x6p linh kién: Quan sat menu bar, chon chUtc nang thich hop va sap 
x€p cdc linh kién trén board mach. O buéc nay, néu can thiét ta c6 thé 
dua vao thém cac linh kién mdi, chang han cac 16 (hold) dé bat 6c cho 
board mach. 

**EEKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 7.8) 
Dat kich thuc céc duOng mach va khai bdo s6 16p: 


e Chon Options Global Spacing ... dé dat khoang cach gilfa Track, 
Via va Pad: 


**EKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 7.9) 
e Chon Options Jumper Settings ..., OK hoac View SpreadSheet 
Layer dé bao s6 16p vé. Trong bai thi nghiém nay, ta vé board mach 
2 16p, nén chi c6 TOP va BOTTOM 1a Routing. 
**EEKSORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 
(Hinh 7.10) 


¢ Chon View SpreadSheet Nets dé dat kich thuGc cdc duOng mach: 


*** SORRY, THIS MEDIA TYPE IS NOT SUPPORTED.*** 


(Hinh 7.11) 


Vé mach: OrCAD h6 tro nhiéu nhiéu chién luc vé mach. Dé khai 
thdc t6t cdc chién luc nay doi hdi phai cd mot s6 kinh nghiém nhat 
dinh. Thong thuOng, ngu0i ta hay vé cdc duOng Neudn — Mass vGi cac tU 
loc nguOn thich hgp tru6c, sau dé thuc hién vé tu d6ng tUing phan va 
cudi cung két ndi cdc phan vUra vé. Nén nhG rang, ché d6 vé tu déng chi 
gitip giam th0i gian, khong cé phan mém nao vé mach tu d6ng tét hon 
tu duy cUa con ngU6i. 


(Hinh 7.12) 


MOt sO thong sO can luu y (Ap dung cho cac sinh vién vé va thué gia céng 
mach Luan van tt nghiép tai mOt s6 co sO 6 Tp. HCM): 


¢ DuOng kinh t6i thi€u cUa cdc 16 chan cam va cac via xuyén mach: 
55 mils (1000 mils=1 inches). Cac 16 chan linh kién nén dat thanh 
hinh oval. 

¢ DO rOng ti thi€u cUa cac dUOng mach (track width): 15 mils (hién 
nay tai Tp. HCM m6t s6 céng ty cé thé thUc hién mach vGi d6 rong 
duOng 6 mils). 

e Khoang cach Track to Track nén >15 mils 

© Dé cé thé vé xuyén qua gitta 2 chan IC, phai dat khodang cach Track 
to Pad thich hOp, nén12 15 mils. 


TV chOn 
Tu thiét ké va m6 phong m6t mach dién tl’ bang Multsim 


Vé mach in cho mach vUa thiét ké bang OrCAD. 


